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We report new results of the NNNLO correction to the S-wave quarkonium wave- 
functions at the origin, which also provide an estimate of the resonance cross section 
in ti threshold production at the ILC. 

1 Introduction 

Top quark pair production near threshold will be an important process at the ILC to deter- 
mine the top quark mass m t , decay width T t and the QCD coupling constant a s . Because of 
high precision required for these quantities, the theoretical uncertainty of the cross section 
should be reduced below a few percent level. For this purpose, the NNNLO QCD calculation 
of the cross section is mandatory. 

Recently we computed the NNNLO correction [TJ[2] to the quarkonium wave-functions at 
the origin, which governs the magnitude of the threshold cross section. In this proceedings 
we present an analysis of the combined result of the papers [TJ [2] . For the details of the 
calculation we refer to the original papers. 

The production cross section of a heavy quark pair QQ is related to the two-point function 
of the vector current 7 71 in QCD: 

(?V - s M V) n( g 2 ) = i J d d xe^(n\Tj^x)f(o)\n), (i) 

where j 71 = Qj^Q, q^ = (2m + E, 0) in the center of mass frame of the QQ, and d = 4 — 2e. 
Near the QQ threshold, the two-point function exhibits the bound-state contribution 

u{q) = 2^— (sT7o) +non - pole ' (2) 

where E n is the energy of the bound state with the principal quantum number n and i 
specifies the physical sheet in the analytic continuation. E n and Z n control the position and 
the height of the resonances in the threshold cross section, respectively. 

The heavy quark threshold dynamics is non-rclativistic (NR), so we utilize an effective 
field theory non-relativistic QCD (NRQCD) for the quark (ijj) and anti-quark (x). In 

*Talk given by Y. Kiyo. Preprint numbers ALBERTA-TH Y- 17-07, PITHA07/14, SFB/CPP-07-68, 
TTP07-29. 
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NRQCD the vector current is mapped onto 



J 



6m 2 r ^ 



(3) 



where c v , d v are matching coefficients, having pcrturbative series expansions in a s . Thus 
the two-point function reduces to the one in NRQCD, whose bound-state contribution is 
expressed by the quarkonium wave- function at the origin, ip n (0), 



i I d d xe tEt (n\T[^a l x}(x)[x f ^mm) ^ 2N c (d - 1) 



E n - (E + i 0) 



+ non-pole. (4) 



The pre-factor 2N c (d — 1) is due to spin(g>color<g)spacc degrees of freedom. The relation 
between the residues of the QCD and NRQCD two-point functions is given by 



m \ 3 



X |Vn(0)| : 



(5) 



where the D 2 term in cq.j3]) was replaced by —mE using the equations of motion of the 
NRQCD fields. The wave-function as well as the matching coefficients possess scale depen- 
dence because of their UV and IR divergences characteristic to effective theory calculations, 
which we treat according to the threshold expansion [3]. The physical quantity measured 
in experiments is Z n , a scale- invariant combination of the matching coefficients and the NR 
wave-function. In the next section we present semi-analytical formulae for all the building 
blocks needed to get Z\, and discuss the importance of the NNNLO correction for stabi- 
lizing the pcrturbative result for the quarkonium wave-functions at the origin against scale 
variation. 



2 NNNLO corrections to the wave-function at the origin 

The wave-function at the origin to NNNLO consists of the Coulomb contribution, the non- 
Coulomb potential contribution, and the ultra-soft correction in NRQCD. The Coulomb 
contribution is finite and calculated analytically in [Hd]. The non-Coulomb [T] and ultra- 
soft 21 computations require regularization and renormalization prescriptions, so that they 
are scheme-dependent quantities. We computed them with conventional dimensional reg- 
ularization and divergences are renormalized in MS scheme. Combining all corrections we 
obtain the following numerical formula for the ground-state wave-function: 



l^i(0)| 2 

l4 0) (o)| 2 



1 + OL s {p) 



(5.25 - 0.32n/)L + 0.21 - 0.13 n/ 



a 2 (ji) 



'18.39 



-2.23 n f + 0.07 n 2 f )L 2 + (1.33 - 0.35 n f + 0.02 n))L + 22.60 - 1.23n/ + 0.02n 2 f 



(53.7 - 9.8 ?i/ + 0.6 n) - 0.01 nf) L 3 + ( - 6.7 + 0.6 n f - 0.07 n 2 f + 0.002 nf) L 2 
-(236.6 - 23.9 n f + 0.8n 2 f - 0.01 n) + 15.0 l m ) L - 22.3 L us + 3.0 l m - 1.5 l 2 n 



-21.0 + 5.0n/ - 0.3 nj + 0.004?^ + 0.0015a 3 + — 



(6) 
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where L = In (/x/(mCi?a s (/i))), Ljjs = ln(e 5 / 6 p/(2m«j(ft))), i TO = ln(/i/m), n/ is the 
number of light quark flavors, ajf] is the constant part of the three loop QCD potential, and 
S e is a contribution from the 0(e) terms of the non-Coulomb potentials given by 



12 



6 2 



(7) 



The effect of <5 e is estimated to be an order of magnitude smaller compared to other constant 
terms [I], so we neglect it in our phcnomcnological analysis. The In 2 a s and lna s [El IS] 
logarithmic terms in eq. © have already been known. 

From the divergent part of the wave-function calculation, the corresponding scale depen- 
dence of c 3 is extracted^ The matching coefficient c v reads 



c v = 1 - — a s (m) + 

07T 



35 /j 2 lln/ 125 C(3) 14 In 2 89 511 
27 U m 2 + 27tt 2 9^ 9 54tt 2 ~ 324 



a s (m) 2 



( 43 35/1/ 
V36^ ~ 162tt 



In 



tfj," / 1399 n/ 2818 851n2\ 
' + V 1944tt ~ 405tt 9tt 



5co 



a s (m) 3 . (8) 



The constant part, <5c3, is not fully known up to now, but the fermionic correction was 
calculated in ITOl, 



t>C 3j n{ = 71/ Cf I> 



39.6CU + 46.7 C F -n f T F — + — - T F — 



163 4tt' 



162 27 



557 26tt' 



162 81 



(9) 



The coefficient d v is known from |llj . and given by 



d v = 1 



16 / u 2 
— l + 31n-^ 
97T V m z 



(10) 



3 Residue of the QCD two-point function 

Now we combine all pieces and show numerical formulae for the residue of the QCD two- 
point function. We use the same coupling a s (/x) for the matching coefficient and the 
NRQCD wave- function to construct the scale- invariant physical residue Z n . 



a Only a Pade estimate 1121 p a[ j e = 6240 (forn/ = 4), 3840 (for rif = 5) is known. 

b The result of [8] has been checked and one term (+ typos) of C3 was corrected in [2]. 

c In cq.JH} a 3 (m) is re-expressed by c* s (/i) using a s (rn)/a a (^t) = 1 + (3q In j^j + 

( "ll^ ) 2 ($0 ln2 ^ + 01 ln ^3") " l where ft are the coefficients of the QCD /3-function in MS-scheme, 

(n* =4,5) 

and ct a = ct s for the bottom and top quarks, respectively. 
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Residue Z,(tt) Residue Z, (bb) 




u [GeV] „ [GeV] 

Figure 1: The scale dependence of the residue of the two-point function for the toponium 
(left) and bottomonium (right), normalized by its zeroth order value at // = mCj?a s (/j). 
The lines refer to LO (black dotted), NLO (red dashed), NNLO (green dashed) and the 
NNNLO (blue solid) for the toponium and bottomonium. 



For the ground state of top and bottom quarkonia, the residue is given by 



ZiS{tt) 



3.66L- 2.13 



8.93 L 2 - 6.14 L + 10.46 - 7.26 l r 



a 2 s (fx) 



18.17 L 3 - 20.26 L 2 + (110.82 - 11.57 l m ) L - 22.27 L us - 16.35 l 2 m - 22.65 l Tl 



(22.60 + 0.0015 a 3 + 0.32 8 e + 0.0645 Sc 3 ) 



(11) 



S(bb) 



3.98 L - 2.00 



a s (n) 



10.55 L - 6.51 L + 11.19 - 7.44 / 



23.33 L 3 - 23.12 L 2 + (125.14 - 14.59 l m ) L - 22.27 L us - 17.36 l 2 m - 26.61 l r , 



(17.44 + 0.0015 a 3 + 0.32 S t + 0.0645 5c 3 ) 



(12) 



where l^isfQQ) (^) | 2 = ( m CF«i5(M)) 3 /(87r) is the LO Coulomb wave-function. To see the 
numerical significance we substitute the following values in the formulae: for the top quark, 
n%t = 175 GeV, // = mt Cf ck s (m) — 32.62 GeV; for the bottom quark, mb = 5 GeV, // = 
nib Cf Oi s {^i) = 2.02 GeV. We use a 3 = a3,p a d , and the unknown 0(e) potentials as well as 
non-n/ term of 5c 3 are set to zero. We obtain the following numbers for the toponium and 
bottomonium ground state at /i = mCFa s (n), 



Zis(tt) 



Z 



(C F mt a s ) 3 



(C F m b a s 



1 - 2.13 a s + 22.7 at + - 38.8 + 5.8 a3 + 37.6 c3 „ ; ct 



lS{bb) 



1 -2.00 a s + 17.9 a z s + - 8.8 + 9.4 a3 + 30.3 c3 ,„i < 



, (13) 
(14) 



where the coupling constant is a s — 0.14, 0.304 for the top and bottom quarkonia, respec- 
tively. 
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In FigfT] we show the scale dependence of the ground-state pole residue for toponium 
and bottomonium. For the NNNLO lines 803 is set to zero, while the gray band indicates 
the size of the contribution from the constant part of C3; the upper/lower edge of the band 
is obtained by taking fermionic corrections £03, n ,/ — 603, nf as an estimate of 8C3 E We 
observe that the scale dependence of the toponium wave-function is reduced significantly at 
NNNLO compared to NNLO as was also observed in renormalization group improved NNLO 
calculations [T31 E]. Its precise value will be fixed only once the third-order matching 
coefficient is completely known. Since the threshold cross section is dominated by the 
ground-state contribution, we expect that the scale dependence of the ti threshold cross 
section will be also improved at NNNLO. For the bottomonium wave- function, strong scale 
dependence remains even at NNNLO and the perturbative expansion may be out of control. 
Only if the constant part of the matching coefficient 803 is negative in total, the scale 
dependence of the bottomonium wave-function at the origin might be acceptable. The 
complete knowledge of C3 is thus mandatory to draw the final conclusion on the size of 
NNNLO correction. 
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d By looking at constant part of the c„ , the non-fermionic correction is larger than the fermionic correc- 
tion in magnitude and the sign is opposite. With this observation, a naive guess for C3 is that the NNNLO 
line in the figure is most likely to be shifted down when the full constant part of C3 is taken into account. 
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